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Abstract 

Using a unified macroscopic QED formalism, we derive an integral equation for the 
van der Waals energy of a two-level atomic system near a carbon nanotube. The 
equation is valid for both strong and weak atom-vacuum-field coupling. By solving 
it numerically we demonstrate the inapplicability of weak-coupling-based van der 
Waals interaction models in a close vicinity of the nanotube surface. 
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1 Introduction 



Carbon nanotubes (CNs) are graphene sheets rolled-up into cylinders of ap- 
proximately one nanometer in diameter. Extensive work carried out world- 
wide in recent years has revealed the intriguing physical properties of these 
novel molecular scale wires [1]. Nanotubes have been shown to be useful for 
miniaturized electronic, mechanical, electromechanical, chemical and scanning 
probe devices and materials for macroscopic composites [2]. Important is that 
their intrinsic properties may be substantially modified in a controllable way 
by doping with extrinsic impurity atoms, molecules and compounds [3]. Re- 
cent successful experiments on encapsulation of single atoms into single-wall 
CNs [4] and their intercalation into single-wall CN bundles [3,5] as well as 
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the theory of spontaneous decay in such systems [6,7] stimulate an in-depth 
analysis of atom-CN van der Waals (vdW) interactions. 

The relative density of photonic states (DOS) near a CN effectively increases 
due to the presence of additional surface photonic states coupled with CN elec- 
tronic quasiparticle excitations [6]. This causes an atom- vacuum- field coupling 
constant (which is proportional to the photonic DOS) to be very sensitive to 
an atom-CN-surface distance. In particular [7], when the atom is close enough 
to the CN surface and the atomic transition frequency is in the vicinity of the 
resonance of the photonic DOS, the system shows strong atom-vacuum-field 
coupling giving rise to rearrangement (" dressing" ) of atomic levels by vacuum- 
field interaction. If the atom moves away from the CN surface, the atom-field 
coupling strength decreases, smoothly approaching the weak coupling regime 
at large atom-surface distances since the role of surface photonic states dimin- 
ishes causing the relative density of photonic states to decrease. This suggests 
strictly nonlinear atom-field coupling and a primary role of the distance- 
dependent (local) photonic DOS in the vicinity of the CN, so that vacuum- 
QED-based vdW interaction models as well as those based upon the linear 
response theory (see [8] for a review) are in general inapplicable for an atom 
in a close vicinity of a carbon nanotube. 

To give this issue a proper consideration, we have developed a simple quantum 
mechanical approach to the vdW energy calculation of a (two-level) atomic 
system in the vicinity of an infinitely long single-wall CN. The approach is 
based upon the perturbation theory for degenerate atomic levels (see, e.g., [9]), 
thereby allowing one to account for both strong and weak atom-vacuum-field 
coupling regimes. In describing the atom-field interaction, we follow a unified 
macroscopic QED formalism developed for dispersing and absorbing media in 
Ref. [10] and adapted for an atom near a CN in Refs. [6,7]. The formalism 
generalizes the standard macroscopic QED normal-mode technique by repre- 
senting mode expansion in terms of a Green tensor of the (operator) Maxwell 
equations in which material dispersion and absorbtion are automatically in- 
cluded. In more detail, the Fourier-images of electric and magnetic fields are 
considered as quantum mechanical observables of corresponding electric and 
magnetic field operators. The latter ones satisfy the Fourier-domain operator 
Maxwell equations modified by the presence of a so-called operator noise cur- 
rent written in terms of a bosonic field operator and an imaginary dielectric 
permittivity. This operator is responsible for correct commutation relations of 
the electric and magnetic field operators in the presence of medium-induced 
absorbtion. The electric and magnetic field operators are then expressed in 
terms of a continuum set of operator bosonic fields by means of the convolu- 
tion of the operator noise current with the electromagnetic field Green tensor of 
the system. The bosonic field operators create and annihilate single-quantum 
electromagnetic medium excitations. They are defined by their commutation 
relations and play the role of the fundamental dynamical variables in terms 
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of which the Hamiltonian of the composed system "electromagnetic field + 
dissipative medium" is written in a standard secondly quantized form. 

Using the approach summarized above, we derive an integral equation for the 
vdW energy of a two-level atomic system near a CN. The equation is repre- 
sented in terms of the local photonic DOS and is valid for both strong and 
weak atom-field coupling. By solving it numerically, we demonstrate the inap- 
plicability of weak-coupling-based vdW interaction models in a close vicinity 
of the CN. 



2 The van der Waals energy of an atom near a carbon nanotube 



Consider a neutral atomic system with its centre-of-mass positioned at the 
point ta near an infinitely long single-wall CN. Assign the orthonormal cylin- 
dric basis {e r ,e^,e 2 } with e z directed along the CN axis. The total nonrela- 
tivistic Hamiltonian of the whole system can then be represented in the form 
(electric-dipole approximation, Coulomb gauge, CGS units) [7,8] 

H = H F + H A + H% + H A 2) F , (1) 
where 

oo 

H F = J dojhw J dR/t(R, u;)/(R, u), (2) 
o 

g A = E — + -E, m ( 3 ) 

H { 21 = - E — Pi- Mr a) + d- Vy3(r A ), (4) 
rrtiC 

i 1 

"S-E^A'M (5) 



are, respectively, the Hamiltonian of the vacuum electromagnetic field modi- 
fied by the presence of the CN, the Hamiltonian of the atomic subsystem, and 
the Hamiltonian of their interaction (separated into two contributions accord- 
ing to their role in the atom-CN vdW interaction - see below). The operators 
P and / in Eq. (2) are those creating and annihilating single-quantum electro- 
magnetic excitations of bosonic type in the CN and the inner integral is taken 
over the CN surface assigned by the vector R = {R cn , <fi, Z} with R cn being 
the CN radius. In Eqs. (3)-(5), rrii, qi, r, and are, respectively, the masses, 
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charges, canonically conjugated coordinates (relative to ta) and momenta of 
the particles constituting the atomic subsystem, d = J2i <?A is its electric 
dipole moment operator. The vector potential A and the scalar potential ip of 
the CN-modified electromagnetic field are given for an arbitrary r = {r, ip, z] 
in the Schrodinger picture by 

oo oo 

A(r) = J duu(iujy l i±(r,uj) + h.c, -V0(r) = J du E"(r, u) + h.c, (6) 





where E (r, uS) = J dr' d ±( -^(r — r') • E(r', a;) is the transverse (longitudinal) 
electric field with ^(r) = -V Q V^(47rr)- 1 and 8^(r) = 5 aP 5(r) - ^(r) 

being the longitudinal and transverse dyadic (^-functions, respectively, and E 
representing the total electric field operator which satisfies the following set 
of Fourier-domain Maxwell equations 

VxE(r,w) = jfcH(r,w), (7) 

- - 47T- 

V x H(r,w) = -ifcE(r,w) H I(r,cA 

c 



Here, H stands for the magnetic field operator, k — cu/c, and 

I(r,u) = JdR5(r-R)J(R,u) = 2J(R cn ,<p,z,u)5(r-R cn ), (8) 



where 



J(R, W )=W /(R, e 2 (9) 



is the operator noise current density associated with CN-induced absorp- 
tion [7], a zz is the CN surface axial conductivity per unit length (in describing 
CN electronic properties, we use the surface axial conductivity model, neglect- 
ing the CN azimuthal current and radial polarizability [11]). 

From Eqs. (7)- (9) it follows that 
47r r 

E(r, u) = i—k dK G(r, R, u) • J(R, u) (10) 



(and H = {ik)~ l V x E accordingly), where G is the Green tensor of the 
classical electromagnetic field in the vicinity of the CN. The set of Eqs. (1)-(10) 
forms a closed electromagnetic field quantization formalism in the presence 
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of dispersing and absorbing media which meets all the basic requirements of 
a standard quantum electrodynamics [10]. All information about medium (the 
CN in our case) is contained in the Green tensor G whose components satisfy 
the equation 

EfVxVx- k 2 )^ G az (r, R, oo) = 8{v - R), (11) 

a=r,tp,z 



together with the radiation conditions at infinity and the boundary conditions 
on the CN surface. This tensor was derived and analysed in Ref. [7]. 

Starting from Eqs. (l)-(5) and using Eqs. (6), (9), (10), one obtains the follow- 
ing secondly quantized Hamiltonian for a two-level atomic system interacting 
with the CN-modified electromagnetic field 

H = Jdoohoo J dR f\R, oo)f(R, oo) + ^a z (12) 



oo 



+ jdoojdR [g (+) {r A , R, oo) ^ - g { ~\r A , R, oo) a] f(R, oo) + h.c. 
o 



Here, the Pauli operators a z = \u)(u\ — a — \l)(u\, = \u)(l\ describe 

electric dipole transitions between the two atomic states, upper \u) and lower 
|/), separated by the transition frequency oo A - This (bare) frequency is modified 
by the interaction (5) which, being independent of the atomic dipole moment, 
does not contribute to mixing the \u) and \l) states, giving rise, however, to 
the new renormalized transition frequency 

ua=u a M -jj^i ]doojdR\g L (v A ,R,oo)A (13) 



in the second term of Eq. (12). On the contrary, the interaction (4), being 
dipole moment dependent, mixes the \u) and \l) states, yielding the third 
term of the Hamiltonian (12) with the interaction matrix elements 

g (±) (r A , R, oo) = g ± (r A , R, oo) ± — g" (r A , R, oo), (14) 



where 

g^\r A ,R,oo) = -i^A^nfkoRea zz (co) d a G^\r A , R, co) (15) 
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with G^(r A ,R,w) = Jdr5^ lu (r A - r) G^ ll) (r, R,oj). The matrix element 
(15), being squared and integrated over the nanotube surface, may be repre- 
sented in the form 

/rfR|g^lD(r A ,R,c)| 2 = (^) 2 ^ID(r A ,c), (16) 

where ^ _L< -"- ) is the transverse (longitudinal) local photonic DOS defined by [7] 

(17) 



with T ±( -^(rA,uj) = Md a dp\m L ^Ga$\YA,YA,u)/hc 2 being the transverse 
(longitudinal) atomic spontaneous decay rate near the CN, and To repre- 
senting the same quantity in vacuum where ImG^a = Sapuj/Qnc. Note that 

Im _L( -"- ) G^"' > = \mG Q a p + Inr L( -"- ) G^jj^ with the second term representing the 
"pure" CN contribution to the total imaginary Green tensor, so that Eq. (17) 

may also be written in the form ^- L (^(r j4 , u) = 1 +£ L ^\ta.,u). 

The Hamiltonian (12) is a starting point for the vdW energy calculation. The 
latter is nothing but the r^-dependent contribution to the ground-state energy 
given by this Hamiltonian. Important is that the two-level atomic subsystem 
(the second term in Eq. (12)) is now described by the renormalized frequency 
(13) which is seen (see Eq. (16)) to decrease with increasing (i.e. when the 
atom approaches the CN surface [6,7]), thereby bringing the two atomic levels 
together, or even making them degenerated if is large enough. To take this 
fact into account in a correct way, one has to calculate the energy using the 
perturbation theory for degenerated levels (see, e.g., [9]). In so doing, one also 
has to account for the upper state degeneracy of the whole system with respect 
to R and uo. In view of this, the wave function of the whole system should be 
written in the form 

oo 

W = C I |OI{0}> + JduJdRC u (R,uj)\u)\{l(R,uj)}), (18) 
o 



where |{0}) is the vacuum state of the field subsystem, |{1(R, oo)}) is its 
single-quantum excited state, Ci and C u are unknown mixing coefficients of 
the lower and upper states of the whole system. The total ground-state energy 
E is then given by the solution of a secular equation obtained by applying the 
Hamiltonian (12) to the wave function (18). This yields the integral equation 



oo 



fi^A f, f Jn \g ( >(r A ,R, uj) 



E = ~^- jdw jdR ^ v ' " , (19) 



o Jtw + — E 
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which the vdW energy E vw is determined from by means of the obvious relation 
E = —fuv A /2 + E vw (r A ) where the first term is the unperturbed ground- 
state energy of the two-level atom. Using further the dimensionless variables 
e vw = E vw /2 r y and x = ftw /2^ (70 = 2.7 eV is the carbon nearest neighbor 
hopping integral appearing in a zz in Eq. (9)), one finally arrives, in view of 
Eqs. (13)-(16), at the following integral equation 



h 



dxT (x)£ (r A ,x) 



(20) 



h 



37T7Q 



r (x) 



dx- 



X + xa - 



h 



J dxr (x)^ ± (r A ,x) - e vw (r A ) 



describing the vdW energy in terms of the distance-dependent photonic DOS 
analyzed in great detail in Ref. [7] and shown in Figure 1(a) for a particular 
case of the atom outside the (9,0) CN. Eq. (20) is universal in the sense 
that it covers both strong and weak atom-field coupling regimes. These two 
limiting cases are obtained by putting x A ~ (V^To^a) Io°dxT (x)^ (r A ,x) 
and e vw ~ 0, respectively, in the denominator of the second term. The former 
is nothing but a degeneracy condition of the unperturbed atomic levels due 
to the atom-field interaction, taking place at small atom-CN-surface distances 
when the local photonic DOS is large enough. The latter occurs at lower 
photonic DOS values when the atom is not very close to the CN surface. In this 
latter case, Eq. (20) gives rise to an ordinary second order perturbation theory 
result with the total Hamiltonian (12) where the last term is considered as a 
perturbation. This result is slightly " improved" compared with those obtained 
from various weak-coupling-based models discussed in the literature (see [8] 
and Refs. therein) since a part of the atom-field interaction is included in the 
unperturbed Hamiltonian via the renormalized transition frequency. 



3 Numerical results and discussion 



Using Eq. (20), we have simulated the vdW energy of the atom outside the 

(9,0) CN. The local photonic DOS functions ^'"(r^x) = ^^(r^, x) — 1 were 
computed in the same manner as it was done in Ref. [7]. The free-space spon- 
taneous decay rate was approximated by the expression T (x) ~ a 3 2'~f x/h 
(a = 1/137 is the fine-structure constant) valid for hydrogen-like atoms [9]. 
Figure 1(a) shows the transverse local photonic DOS, as a typical example, for 
the atom located at several distances from the CN surface. The DOS function 
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Fig. 1. Typical example of photonic DOS (the transverse local photonic DOS) (a) 
and the vdW energy (b) of the atom outside the (9,0) CN. The (bare) atomic 
transition frequencies are indicated by dashed lines in Fig. 1(a). 

is seen to increase with decreasing the atom- CN- surf ace distance, representing 
the increase of the atom-field coupling strength as the atom approaches the 
nanotube surface [7]. The vertical dashed lines indicate the bare atomic tran- 
sitions frequencies xa for which the vdW energies shown in Figure 1(b) were 
calculated. The frequencies are typical for heavy hydrogen-like atoms such as 
Cs for which xa may be estimated from its first ionization potential [12] to be 
3.89 eV x 3/4 x (270) _1 ~ 0.5 (the factor 3/4 comes from the Lyman series 
of H), or less for highly excited Rydberg states. 

In Figure 1(b), the vdW energies given by exact numerical solutions of Eq. (20) 
are compared with those obtained from the same equation within the weak 
and strong coupling approximations. At small atom-CN-surface distances, the 
exact solutions are seen to be fairly well reproduced by those obtained in the 
strong coupling approximation, clearly indicating the strong atom-vacuum- 
field coupling regime in a close vicinity of the nanotube surface. The devia- 
tion from the strong coupling approximation increases with xa, that is easily 
explicable since the degeneracy condition of the unperturbed atomic levels is 
more difficult to reach for larger interlevel separations. As the atom moves 
away from the CN surface, the exact solutions deviate from the strong cou- 
pling solutions and approach those given by the weak coupling approximation, 
indicating the reduction of the atom-field coupling strength with raising the 
atom-surface distance. The weak-coupling solutions are seen to be divergent 
close to the nanotube surface as it should be because of the degeneracy of the 
unperturbed (bare) atomic levels in this region. 

To conclude, we have investigated vdW interactions of a two-level atomic sys- 
tem near a single-wall CN within the unified QED approach accounting for 
both strong and weak atom- vacuum-field coupling. We have demonstrated the 
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inapplicability of weak-coupling-based vdW interaction models in a close vicin- 
ity of the nanotube surface where the local photonic DOS effectively increases 
giving rise to atom-field coupling enhancement followed by the degeneracy and 
rearrangement ("dressing") of bare atomic levels. 

The work was performed within the framework of the Belgian PAI-P5/01 
project. LB. thanks the Belgian OSTC. 
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